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Abstract 

The propagation of longitudinal acoustic waves in multilayer structures based on porous silicon and the experimental 
measurement of acoustic transmission for the structures in the gigahertz range are reported and studied theoretically. 
The considered structures exhibit band gaps in the transmission spectrum and these are localized modes inside the 
band gap, coming from defect layers introduced in periodic systems. The frequency at which the acoustic resonances 
appear can be tuned by changing the porosity and/or thickness of the defect layer. 
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Background 

The study of acoustic and elastic wave propagation in 
phononic crystals (PCs) [1-3] have been studied theo- 
retically [4] and experimentally [5] in recent years. In 
analogy with the photonic band gap materials, empha- 
sis in phononic crystals has been on achieving large 
acoustic band gaps within which propagation of sound is 
forbidden. 

The interest about phononic band gap materials is in 
both, pure and applied physics. Fundamental physics has 
a special interest concerned with the localization phe- 
nomena of sound and vibrations in PCs. Researchers have 
prospected numerous applications based on cavity struc- 
tures built around PCs, such as wave filters, waveguides, 
and splitters [6-9]. Furthermore, it is possible to design 
cavities for coherent (single-wavelength) phonon gener- 
ation and control, to attain phonon amplification and 
'lasing' in the called 'saser', one of the most important 
potential applications [10-12]. 

Periodic solid-state structures exhibit transmission stop 
bands for waves at certain frequencies. By placing one or 
more defects into a perfect phononic crystal, acoustic cav- 
ities are created inside the system. The presence of these 
defects, produces localization of elastic or acoustic modes 
inside the phononic band gap. These localized modes are 
the acoustic analog of donor or acceptor states produced 
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inside the band gap of semiconductors. In analogy with 
electronic systems, one can consider these acoustic states 
to control the sound propagation through the structure. 

If a defect is introduced into a periodic structure, the 
translational symmetry is broken and highly localized 
defect modes within the band gaps are created [6,8,13,14]. 
Point, linear, and planar defect states have been the- 
oretically investigated in one-dimensional (ID), two- 
dimensional (2D), and three-dimensional (3D) phononic 
crystals [3,15,16]. In ID structures, a microcavity can be a 
spacer layer of thickness A. /2 enclosed by two Bragg reflec- 
tors [17]. In 2002, Trigo et al. proposed phonon cavities 
in structures consisting of two semiconductor superlat- 
tices enclosing a spacer layer, showing that acoustical 
phonons can be confined in such layered structures if the 
spacer thickness is an integer multiple of the acoustic half- 
wavelength at the center of one of the superlattice-folded 
minigaps. These acoustic cavities are semiconductor mul- 
tilayers in the nanometer scale and are fabricated by 
molecular beam epitaxy (MBE), which is a sophisticated 
and expensive technique that requires ultra-high vacuum 
system and a very tight control on the growth parameters, 
and modulate the thicknesses is easier than to modulate 
the elastic properties of the layers. Contrasting, porous 
silicon (PS) multilayer fabrication is relatively easy and 
considerably less costly, besides that this material allows to 
modulate both the thicknesses and the elastic properties 
of each layer. 

PS is known as a versatile material with applications in 
light emission, sensing, and photonic devices [18]. The 
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possibility of producing acoustic band gaps in PS was 
proposed in 2005 [19], and detailed calculations of pre- 
dicted bandwidths were subsequently published [20]. 
Recently, experimental results of Brillouin light scatter- 
ing suggested the existence of zone-folded phonons and 
phononic band gaps in PS multilayers [21]. G. N. Aliev 
et al. in 2010 reported experimental measurements of 
transmission for longitudinal acoustic waves and observa- 
tion of hypersonic band gaps in a periodic PS structure 
with an acoustic microcavity consisting of a spacer layer 
of thickness k/2 [22]. 

By taking advantage of the possibility to modulate the 
elastic properties of PS layers, and considering that it is 
possible to create localized modes by introducing a defect 
layer with different acoustic properties into a periodic 
structure, in this paper, we investigate the propagation of 
longitudinal acoustic waves in multilayer structures based 
on PS, that exhibit resonant cavity modes in frequencies of 
gigahertz (GHz), consisting of defect layers intentionally 
introduced in periodic structures. The design and mate- 
rial parameters that allow to create these localized acous- 
tic modes is discussed, and experimental results of the 
measured acoustic transmission in PS samples fabricated 
by electrochemical etching are presented. 

Methods 

Theoretical models 

The multilayer PS structures studied here have thick- 
nesses in micrometer range and the procedure used to 
fabricate them creates mesoporous silicon with an aver- 
age pore diameter of 20 to 50 nm. On the other hand, 
in our experiments, the typical longitudinal wavelengths 
excited throughout the samples are 3 to 7 urn depending 
on porosity. Accordingly, each of the individual layers in 
the structures is assumed to be homogeneous. 

The longitudinal acoustic wave equation in the contin- 
uum limit for a solid inhomogeneous along the z direction 
(but homogeneous along the x and y directions) is given 
by [23], 



d 2 u(z, t) 



= (k + 2n)j 



d 2 u(z, t) 
dz 2 ' 



(1) 



where pj is the mass density, (A, + 2/x)y = pv 2 j and u(z, t) 
is the atomic displacement. Here, / is an index identifying 
each layer. The limits of the elastic continuum descrip- 
tion of wave propagation in ordered media depends on 
the dimensions of the system compared with the wave- 
length. When the dimensions approach nanometer-length 
scales, atomistic treatments using first principles or semi- 
empirical methods may become necessary [24]. How- 
ever, in our case, the thicknesses of the layers are in the 
micrometer range and each layer can be considered as a 



homogeneous layer; thus, the model described before is 
assumed valid. 

In a solid, the acoustic waves can be longitudinal or 
transversal. In this letter, only longitudinal waves propa- 
gating through PS are considered because in our exper- 
iments, the waves are coupled to the samples through a 
liquid at normal incidence. The mass density p is a func- 
tion of the porosity and is described by p = po{l — P) 
where po = 2.330 g/cm 3 is the density of bulk silicon 
and P the porosity. The acoustic velocity dependence on 
porosity is given empirically by \>i = vlo(1 — P) k > being 
Vlo the longitudinal velocity of sound in bulk silicon along 
the (100) crystallographic direction and k > 0.5 is a con- 
stant [25-28]. In general, the parameter k depends on PS 
morphology which in turn depends on the doping level of 
the Si substrate [25,26]. 

The general solution of Equation 1 takes the form, 

u(Zj, t) = ^4 ; + exp (ikjzj) + AJ exp (— ikjzj) j exp(icot), 

(2) 

corresponding to waves propagating to the right (A + ) 
and left {A~). At each interface, this solution must satisfy 
the boundary conditions related to the continuity of the 
atomic displacement and stress, 



uj {dj) = Uj+\ (dj) 



and 



3 Mi (dj) 



2 9m /'+1 (4) 
Z=dj - Pj+\v L{j+1) — 



(3) 



(4) 



z=dj 



respectively. Here, dj denotes the position of the ;'-th inter- 
face between / and j+1 layers. The frequency a> is related 
to its wave vector via a> = kjVj, with Vj the sound speed 
in the ;'-th layer and w = 2nf, being/ is the frequency 
in s _1 . Using the transfer matrix method (TMM) [26], we 
can relate the amplitudes of the fields A~J and AJ in the 
layer / of the system with the amplitudes of the wave in the 
j+1 layer according to 



A 



A j ) ' \ A j+i 



7+1 



(5) 



The transfer matrix 7) appearing in the previous equation 
propagates the amplitudes through a layer with thickness 
dj, mass density pj, and sound longitudinal velocity vq, 
and is given explicitly by, 



cos 



(kjdj) 



—isin(kjdj) 
k AjP) 

-ikjV 2 jPj sin (kjdj) cos (kjdj) 



(6) 
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If we consider a structure formed by N layers, the total 
transfer matrix representing the structure is obtained by 
multiplying, in the appropriate order, a series of AT transfer 
matrices, each one given by a matrix of the type appearing 
in Equation 6. The obtained matrix relates the displace- 
ment vector at the beginning of the structure with that at 
the end, and represents a 2 x 2 set of equations that can 
be fully solved. 

With the above formalism, one can derive the acoustic 
eigenenergies and eigenvectors. The reflectivity and trans- 
mission can also be calculated as the square modulus of 
Aq and Aft, respectively, by imposing the boundary con- 
ditions Aq = 1 and AJj = 0 for a wave traveling from 
right to left. Here 0 and N label the first and last layer of 
the structure, respectively. 

Attenuation can be included by taking the wave vec- 
tor kj complex, such that Kj = kj — a,-, where a; is 
attenuation coefficient. The form of the attenuation coef- 
ficient depends on the physical process causing loss and 
we assume that the Akhiezer model is dominant in a semi- 
conducting material. This gives a = i]co 2 /2pv 3 , where r) is 
the viscosity. However, it is known that introducing acous- 
tic attenuation into the model leads to important effects 
as the shrinking of gaps, only for frequencies higher than 
180 GHz [29]; therefore, no absorptive behavior is consid- 
ered in our model since no important effects are obtained 
if they are included. 

Furthermore, the position and width of the band gap 
are critical parameters for devices that reflect or localize 
the acoustic waves [30]. Band structures of many kinds 
of periodic phononic crystals have been reported [31-33]. 
The most commonly studied acoustic band gaps in ID 
PCs are the Bragg type, appearing at an angular frequency 
co of the order of VL(T)/d (vl(T) refers to the longitudi- 
nal (transverse) elastic wave velocity and d is the lattice 
constant). 

An acoustic Bragg mirror can be made by repeating n 
times a basic block of two materials with different acous- 
tic properties. The mismatch in the acoustic impedance 
between these layers results in waves that are reflected and 
interfere, giving an acoustic band gap around a central fre- 
quency f B . For normal incidence, this frequency is given 
by 

f B = -( Pl d 1 /Z l + p2d2/Z 2 r 1 , w = l,2,3,..., (7) 
m 

being m the order of the stop band, d\ and d 2 are the layer 
thicknesses, and Z\ and Z 2 are the acoustic impedances of 
layers 1 and 2, respectively. The acoustic impedance Z is 
given by pv, with v as the sound velocity and p as the mass 
density. The condition p\d\IZ\ = p 2 d 2 /3Z 2 optimize the 
stop-band width and reflectivity, corresponding in an infi- 
nite stack, to the first minigap at the Brillouin zone center. 
The reflectivity at the center of the stop-band depends on 



the acoustic impedance mismatch between the two mate- 
rials Z 2 jZ\, and for n pairs of layers is given by [17,22], 



Rb 



(Z2/Z1) 



2h 



1 



LCZ2/Z1) 



2h 



(8) 



In [34], the authors considered periodic semiconductor 
structures of GaAs/AlAs to introduce microcavities as 
spacer layers of thickness k/2. However, for a 10-period 
GaAs/AlAs mirror, R B ~ 0.880, while R B ~ 0.996 if 
n = 20. For a PS structure, a porosity variation of 15% 
between the constituent layers of 52% and 67% of poros- 
ity, leads to R B ~ 0.997 for n = 6. Thus, by modulating 
the porosity of the PS structures, very high reflectivity 
values can be achieved. This is an essential condition to 
obtain narrow transmission bands into the stop bands 
corresponding to the cavity modes. 

To demonstrate the localization in time domain, we con- 
sider the propagation of a Gaussian pulse through the 
structure. The Gaussian pulse is described by g(f) = 
exp ^— 47T [(/ —fo) / cr ] 2 ^ were fo is the central frequency 
and a the pulse width. In response to the incident pulse, 
the time and spatial variations of the displacement field 
u(z, t) inside the sample can be calculated according to the 
scattering state method as [35], 



1 r + °° 

u{z,t) = — / u(z,f)g(f)e- l2n fidf 

27T J-oo 



(9) 



where u(z,f) is the displacement field distribution at 
each frequency, which is obtained by the transfer matrix 
method. 

Experimental details 

Samples were electrochemically etched from boron- 
doped (lOO)-oriented Si substrates with a resistivity 
of 0.007 to 0.013 Q.cm. Room-temperature anodization 
was performed using a 1:1 solution of HF (40%) and 
ethanol (99.98%). The acoustic transmission measure- 
ments reported here were done using a Vector Network 
Analyzer (VNA). Each sample was placed between two 
ZnO-based piezoelectric transducers with a central fre- 
quency of 1.1 GHz and an operation bandwidth of 500 
MHz. The transducers consist of a piezoelectric layer driv- 
ing waves into a silicon pillar with a thickness of 520 [im. 
To couple the transducers to the specimen, In-Ga eutec- 
tic was used. The transducer front surface was aligned 
parallel to the sample surface using two orthogonal micro- 
scopes so that the acoustic waves impinge normally into 
the PS layers. The transducers were connected to the 
VNA ports and transmission parameters were measured 
as function of frequency, more details of the experimental 
set-up can be found in [36]. Dependence of sound veloc- 
ity with porosity for Si wafers and experimental conditions 
used here were determined previously and reported in 
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[36], the values of viq and the parameter k used here are 
8.44 km/s and 0.56, respectively. 

For the fabrication of PS multilayers, we consider the 
inclusion of 'etch stops' or 'etch breaks' where the current 
is interrupted to stop the etching of the Si wafer in order 
to prevent any depletion of HF [37]. The introduction of 
these etching breaks is necessary to obtain layers with 
constant porosity with depth [38]. Because our samples 
include very thick layers, with large mismatch porosi- 
ties between them, the number and length of the etch 
breaks are important to obtain homogeneous structures. 
We found that etch breaks of 4 s with a ratio (etch break 
time) /(etching time) from 3.3 for low porosities (52%) to 
7.3 for high porosities (67%) are enough to minimize any 
chirp in the layers. 



Results and discussion 

Thicknesses of layers were measured by optical 
microscopy, and the layer porosities were determined 
from optical reflectance spectra by fitting our experi- 
mental measurements and comparing them with our 
theoretical simulations for each sample. The acoustic 
transmission and field intensity distribution have been 
modeled using the transfer matrix method described 
before and taking into account the effect of the sam- 
ple (PS-Si substrate), transducers (Si pillars), and In-Ga 
eutectic liquid used to couple the transducers to the 
sample. 

Three PS multilayer samples are considered here to 
show the effect of localization inside the structures. All 
of them consist of layers a and b repeating alternatively, 




Substrate 




< 

Depth in the sample (urn) 



Figure 1 Acoustic transmission and distribution of the displacement field for the periodic case, sample 1. (Top) Scheme of the periodic 
structure consisting of 12.5 periods of layers a and fa. (Middle) Acoustic transmission spectra, measured in solid line and calculated in dashed line. 
The measured transmission, recorded on a logarithmic scale, is normalized to its maximum and corrected by an envelope function of the transducer 
response. (Bottom) In solid line, squared phonon displacement corresponding to the central frequency of the gap. The dashed line represents the 
material acoustic impedance and serve to identify the position in the sample. 
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and a defect layer, c, in the middle of the structure. The 
sequence used for structures was abababababab — c — 
babababababa = (ab) 6 c(ba) . In the first sample (1) 
porosities and thicknesses of layers a and b are P a = 53%, 
d a = 1.15 urn, Pb = 67%, and db = 1.10 um, respectively. 
Here, layer c has the same thickness and porosity of layer 
a, and therefore, this sample is completely periodic. The 
porosities and thicknesses of the layers were chosen to 
obtain the fundamental stop band within the bandwidth 
of the acoustic transducers, and satisfying Equation 7. 
A scheme of structure 1 is displayed in the top panel 
of Figure 1. The central panel of Figure 1 (solid line) 
shows the measured acoustic transmission spectrum of 
the PS periodic structure with a total thickness of approxi- 
mately 27 urn. The band gap in the transmission spectrum 



observed around 1.15 GHz and ranged from 0.94 to 1.38 
GHz is the first-order acoustic stop band of the mirror, 
corresponding to m = 1 in Equation 7. This fundamen- 
tal stop band shows an attenuation of approximately 50 dB 
with a fractional bandwidth of 37%. The dashed line is 
the result of calculations using TMM. Good agreement 
between modeled and measured spectra is seen. The fine 
features of the spectrum are not noise but the longitudi- 
nal modes of the Si pillars of the transducers and the Si 
substrate of the sample. The fundamental band gap has a 
depth of approximately 50 dB which is less than the mod- 
eled value of approximately 100 dB. However, this is an 
experimental limitation on the setup used and its noise 
floor. Thus, we can conclude that the stop band has a 
depth of at least 50 dB. The bottom panel of Figure 1 




0.6 0.8 1.0 1.2 1.4 1.6 1.8 
Frequency (GHz) 




Depth in the sample (|um) 

Figure 2 Acoustic transmission and distribution of the displacement field for sample 2. (Top) Scheme of a structure consisting of two mirrors 
with six periods of layers a and fa enclosing a defect layer of higher porosity between them. (Middle) Measured acoustic wave transmission 
spectrum through the sample (solid line). The dashed curve is the calculated spectrum (see text for details). (Bottom) In solid line, squared phonon 
displacement corresponding to the first (thick line) and second (thin line) modes at 1 .01 and 1 .27 GHz, respectively. The dashed line represent the 
layer acoustic impedance. 
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shows the squared displacement field corresponding to 
the central frequency of the gap, 1.15 GHz. The dashed 
line represents the material acoustic impedance and is 
useful to identify the position in the sample. As can be 
seen, the displacement field is not localized, as is expected. 

Now, based on the concepts mentioned before about 
cavities, we will show how the intentional introduction of 
a defect layer between a pair of mirrors can lead to for- 
mation of an acoustic cavity mode within the stop band. 
For this purpose, we consider two structures: sample 2 and 
sample 3. 

In sample 2, porosities and thicknesses of layers a, b, 
and c are: d a = 1.15 um, P a = 52%, dy = 1.00 urn, 
P b = 65%, d c = 1.15 urn and P c = 74%, respec- 
tively. The defect (layer c) corresponds to a layer with the 
same thickness, as the periodic case, but higher poros- 
ity (lower impedance), as is shown schematically at the 



top of Figure 2. In the middle of Figure 2 are shown the 
acoustic transmission spectra, measured experimentally 
(solid line) and calculated theoretically (dashed line). The 
introduction of the defect layer results in well-localized 
transmission modes at 1.01 and 1.27 GHz, within the fun- 
damental stop band ranged from 1.02 to 1.47 GHz and 
with a fractional bandwidth of 35%, as it can be seen in 
the transmission spectrum. At the bottom of the Figure 2 
is shown (in solid line) the displacement field distribution 
as a function of the position in the sample, correspond- 
ing to the cavity modes, the first (thick line) and second 
(thin line) modes at 1.01 and 1.27 GHz, respectively. It can 
be seen that the amplitude of the acoustic displacement is 
maximum around the defect layer. The dashed line is the 
material acoustic impedance. 

Sample 3, represented schematically at the top of 
Figure 3, contains a defect consisting of a layer with lower 
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Figure 3 Acoustic transmission and distribution of the displacement field for sample 3. (Top) Scheme of a structure of two mirrors with six 
periods of layers a and fa enclosing a defect layer of lower porosity. (Middle) Measured (solid line) and calculated acoustic transmission spectra (see 
text for details). (Bottom) In solid line, squared phonon displacement corresponding to the cavity mode frequency (thick line) at 1 .1 5 GHz, and for a 
frequency of 1 .46 GHz (thin line). The dashed line represent the layer acoustic impedance. 
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porosity (higher impedance) at the center of the struc- 
ture. Here, thickness and porosities are: d a = 0.89 um, 
P a = 65.5%, d b = 1.12 um, P h = 53%, d c = 0.89 um, 
P c = 42%, for layers a, b, and c, respectively. The defect 
layer (c) keeps the periodicity in thickness but the porosity 
changes. As it can be clearly seen in measured trans- 
mission spectrum shown in Figure 3, this results in an 
acoustic cavity mode at 1.15 GHz within the fundamental 
stop band ranging from 1.02 to 1.44 GHz (34% fractional 
bandwidth). The corresponding displacement field distri- 
bution for this cavity mode is shown at the bottom of 
the same figure (thick line) and demonstrates that the 
displacement field is maximum around this cavity in the 
same way as the second mode in sample 2. For demonstra- 
tion purposes, we have calculated the displacement field 
for 1.46 GHz and the results are shown in Figure 3 using a 
thin line. Localization effects cannot be observed. 

In Figures 1, 2, and 3, good agreement between mod- 
eled and measured spectra is observed, and the slight 
differences between theoretical and experimental acoustic 
transmissions are due to features of porous silicon layers 



which are not considered here, as the roughness at the 
interfaces, as well as intrinsic error coming from the mea- 
sured procedure, and not to absorption properties, as was 
explained before. 

In Figure 4, we show the time-resolved displacement 
field u(z,t), corresponding to the time evolution of a 
Gaussian pulse in the samples calculated using Equation 9. 
Figure 4a,b corresponds to the time and spatial varia- 
tions of the displacement field inside sample 2, using 
fo = 1.01 GHz in Figure 4a and 1.27 GHz in Figure 4b. 
These values correspond to the frequencies where the 
first and the second cavity modes appear, respectively. 
Figure 4c shows the displacement field inside of sam- 
ple 3 for fo = 1.15 GHz, the frequency of the corre- 
sponding cavity mode. Figure 4d corresponds to sample 
3 using fo = 1.46 GHz. We use a pulse with a = 
200 MHz for all cases. In Figure 4a, it can be seen that 
the displacement field is in the center of the PS struc- 
ture, corresponding to the defect layer. Figure 4b corre- 
sponds to the same structure but the incident Gaussian 
pulse has a frequency equal to 1.27 GHz; in this case, 
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Figure 4 Time evolution of acoustic Gaussian pulses. Calculations of ln(u(z, f)) for acoustic Gaussian pulses centered at the frequency fo 
indicated in each figure and a = 200 MHz for all cases. (a,b) Sample 2. (c,d) Sample 3. Zero in x-axis is placed at the surface of the PS sample. 
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localization is around the defect layer and not inside it. 
In Figure 4c, as Figure 4b, the localization is around the 
defect layer and not inside it, because this corresponds to 
the second mode of sample 3. Calculations for sample 3 for 
the peak at 1.46 GHz appears in Figure 4d, as can be seen, 
localization of the displacement field is not observed. 

In order to estimate the displacement field intensity 
u(z, t) 2 , within the defect layer as a function to the time, 
we integrate the displacement field on the defect layer 
using, 

= J u 2 (z,t)dz, (10) 

J defect 

where $>(t) is the displacement field intensity contained 
in the defect as a function of the time. Figure 5 a shows 
^(t) for sample 2 for two Gaussian pulses centered in the 
frequencies fo indicated there, as expected the first mode 
has higher displacement field intensity in the defect layer 
because the acoustic wave is localized in the center of the 
PS structure, on the contrary to the second mode, see 
Figure 5a. In the case of sample 3, the localization is less 
than sample 2 for the two Gaussian pulses considered, that 
is, the <E> (t) amplitude for sample 2, see Figure 5b, in the 
first mode is around 30 times more the <i> (t) amplitude in 
sample 3 for one incident Gaussian pulse with frequency 
equal to 1.15 GHz. Finally, localization is not observed in 
sample 3 for a Gaussian pulse with a frequency of 1.46 
GHz, as is expected, see Figure 5b. 

The modeled transmittance of the periodic case (sample 
1) and for the two cavity structures (samples 2 and 3), 
obtained by the TMM, shows a good match with the 



experimental results. The localized acoustic resonances 
can be tuned at different frequencies (within the acoustic 
band gap) by changing the porosity of the defect layer. 

Moreover, for commercial acoustic mirrors which are 
components of solidly mounted resonators and filters [39], 
a low-acoustic-impedance material such as Si02 is lay- 
ered with high-impedance materials such as tungsten or 
molybdenum. Following Equation 8, for the layer pair 
of molybdenum and silica, where acoustic impedances 
are 66.2 MRayl and 13.1 MRayl, respectively, the fixed 
impedance ratio is 5.1, and the same impedance ratio can 
be obtained using PS layers of 30% and 75%, so, by mod- 
ulating the porosity, very high reflectivity values can be 
achieved. Besides, by increasing the number of pairs of 
layers enclosing the defect layer leads to narrowing of the 
transmission band together with deepening of the stop 
band. 

Conclusions 

We have demonstrated theoretically by using the TMM 
and experimentally by acoustic transmission measured 
directly, the formation of acoustic cavity modes in GHz 
frequencies by introduction of defects into periodic struc- 
tures based on PS. Acoustic resonances can be tuned 
at different frequencies by changing the porosity of the 
defect. And we proved that these resonant modes appear 
due to the localization of the field into the defect. The 
acoustic mirrors and cavity structures based on PS have 
a performance which is at least comparable with that 
devices based on semiconductor superlattices. This study 
could be useful for the design of acoustic devices, such as 
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Figure 5 Displacement field intensity as a function of time. Theoretical calculations for the displacement field intensity u(z, f) 2 in the defect as a 
function of time for (a) sample 2 and (b) sample 3, for frequencies indicated in each figure. 
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highly selective frequency filters with applications in GHz 
range. 
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